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ABSTRACT 

NGC 4203 is a nearby early-type galaxy surrounded by a very large, low-column-density H i 
disc. In this paper we study the star formation efficiency in the gas disc of NGC 4203 by using 
the UV, deep optical imaging and infrared data. We confirm that the H i disc consists of two 
distinct components: an inner star forming ring with radius from ~1 to ~3 R^//, and an outer 
disc. The outer H i disc is 9 times more massive than the inner H i ring. At the location of the 
inner H i ring we detect spiral-like structure both in the deep g' - r' image and in the 8 yum 
Spitzer-IRAC image, extending in radius up to ~ 3 R^//- These two gas components have a 
different star formation efficiency likely due to the different metallicity and dust content. The 
inner component has a star formation efficiency very similar to the inner regions of late-type 
galaxies. Although the outer component has a very low star formation efficiency, it is similar 
to that of the outer regions of spiral galaxies and dwarfs. We suggest that these differences 
can be explained with different gas origins for the two components such as stellar mass loss 
for the inner H i ring and accretion from the inter galactic medium (IGM) for the outer H i 
disc. The low level star formation efficiency in the outer Hi disc is not enough to change the 
morphology of NGC 4203, making the depletion time of the H i gas much too long. 
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1 INTRODUCTION 


Although we have long known that early-type galaxies (ETGs) can 
harbour a stellar disc, their actual distribution of bulge-to-disc ra¬ 
tios B/D has been debated for many decades (e.g., Spitzer & Baade | 


|1951[ ISandage, Freeman & Stokes]|1970[ [van den Bergh 1976| l. 

Recently, a number of studies based on broad-band imaging and 
integral-field spectroscopy at optical wavelengths have concluded 
that most ETGs host a stellar disc, and that ETGs as a family cover 
a similar range of B/ D ratio as spirals ([Cappellari et al.|2011[ Lau- 


rikainen et al.||2011[ [Kormendy & Bender||201^ Krajnovic et al. 


201 3[ Weijmans et al.|2014| |. Discs are detected also with a num¬ 

ber of additional techniques such as narrow-band imaging of emis¬ 
sion lines (e.g., H/5), imaging of dust absorption, or the observation 
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of the kinematics of planetary nebulae and globular clusters ( |Ro-| 
[manowsky et al.|2()03l|Sarzi et al.|2006||Pota et al.|2013| l. Taken 
together, these results suggest that galaxy evolution is dominated 
by processes which allow the growth and survival of stellar discs in 
galaxies of both early- and late type. 

Discs have also been found in the cold gas phase of ETGs by 
means of both neutral hydrogen (Hi) and CO observations (e.g., 
van Driel & van Woerden|19^|Young|2002[|Morganti et al.|2006[ 

Oosterloo et al.||2007||2010[ [Crocker et al.||201l| . More recently, 

ATLAS^° observations of H i and CO have shown that ~50 percent 
of all ETGs contain some cold gas ( [Young et al.|2011[|Serra et ^ 
|2012| hereafter S12). Molecular gas is typically found in small gas 
discs in the central regions and is linked to small amounts of star 
formation ( jAlatalo et al.||2013[ [Davis et ^[2013| l. In contrast, in 
~20 percent of ETGs outside clusters H i has been found distributed 
in low-column density discs or rings with typical sizes of many 
tens of kpc, much larger than the stellar body (S12). If the column 
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Figure 1. SDSS g-band(left) and new deep WSRT Hi (right: see Sec. HD images of NGC 4203. The size of the images is ~ 74 x 74kpc^ (17' x 17'). The 
colourbar refers to the Hi image. The Hi beam is shown in the bottom-left corner of the Hi image. The blue boxes represent the size of the images shown in 

Fig-El 


density of the gas is high enough, these H i discs are potential fuel 
for star formation (SF). Contrary to the general idea that galaxies 
evolve from late- to early type, these H i reservoirs could generate 
the growth of stellar discs and therefore trigger the transition of the 
host from an early- to a late-type morphology ( |Cortese & Hughes] 
|2009[[Fang et al.|2012||Salim et al.|201^ . 


In this context, many recent studies have focused on the pres¬ 
ence of SF in the outer discs of ETGs ( Salim & Rich|2010'[|Lemo-| 


|nias et al.|2011| [Moffett et al.|2012[ [Salim et al.|2012| ). Such SF is 

typically traced by UV light emitted by massive stars of age up to 
several hundreds Myr ( [Kauffmann et al.|2()07l|Bianchi|201 l[|Salim| 
|et al.|2012] l. |Moffett et al.| ( |2012| ) show that -40 percent of 38 red 
sequence and blue cloud ETGs with stellar mass below ~ 4 x 10^^ 
Mq have extended UV discs, which is evidence for disc growth in 
these objects. 


In this paper, we aim to bring these two lines of investiga¬ 
tion (Hi and UV) together and study SE in the large Hi disc of 
the early-type galaxy NGC 4203. This is a nearby, nearly face-on 
lenticular. Eigure[^hows an SDSS g-band image of the galaxy (left 
panel) and Table provides some basic information. NGC 4203 
has been observed at many wavelengths. Optical studies show that 
NGC 4203 has a relatively massive stellar disc {B/D - 0.4-0.5 - 
|Burstein||1979[ |Krajnovic et al.||2013| ). Early H i observations re¬ 
vealed a low-column-density H i disc with a complex morphology 
( [Burstein & Krumm|1981[|van Driel et al.|1988| l. More recent ob¬ 
servations show that the H i disc extends to very large radius and is 
characterised by spiral-like Hi arms (S12). This gas rotates regu¬ 
larly around the galaxy and is kinematically aligned with the stellar 
disc ( |Serra et al.|2014| ). 

Despite the large Hi reservoir (~ 10^ M©), the stellar popu¬ 
lation is uniformly old (~ 10 Gyr) inside 1 /?eff ( |McDermid et al.| 
|2015| ) suggesting that the H i is not fueling a high level of SE in this 
region. Yet, there are signs of possible low-level SE within the stel¬ 
lar body. Eor instance, NGC 4203 contains 2.5x 10^ M© of H 2 in the 
centre ( [Welch & Sage|2003l|Young et al.|201 l| l. Hubble Space Tele¬ 
scope observations show nuclear dust structures in the very centre 


Table 1. General properties of NGC 4203 


Quantity Value Unit Reference 


Hubble Type 
R.A. (J2000) 

Dec. (J2000) 
Helio. Rad. Vel. 
Distance 

A5 

Effective Radius 
log(Af/ B)star 
logMstar 

Galactic E{B - V) 


SABO 

12 ^ 15 '” 05 . 06 ^ 

-f33'^1U50.382^ 

1086+4 

14.7 

-23.4 

0.05 

29.5 
0.5 

10.6 
0.012 


- 

(1) 

HMS 

(2) 

DMS 

(2) 

km/s 

(1) 

Mpc 

(3) 

mag 

(3) 

mag 

(3) 

arcsec 

(3) 

M©/L©^ 

(4) 

M© 

(4) 

mag 

(5) 


Referenc es:(l) j de Vaucou leurs 


j2005) (3)|Cappellari et al. 
[Finkbeiner & DavisH1998j ;. 


et al.j |T9^; (2) [Anderson & Ulvest^ 

1 ; (4 ]Cappellari et al.[j2013) ; (5) [SchTeg^ 


of the galaxy (300 pc radius; [Erwin & Sparke|2003} . The dust is 
visible also on a large scale in a 8 fim image, where the flux, mostly 
from polycyclic aromatic hydrocarbons (PAHs) appears distributed 
in a spiral-like pattern ( [Pahre et al.|2004| we will make use of this 
image as well as of a 24 fim image in support of our analysis of SE 
in NGC 4203). 

Here we use new, deep H i observations and UV imaging from 
the Galaxy Evolution Explorer (GALEX) to study the relation be¬ 
tween Hi and star formation as well as the impact of this large 
gas reservoir on NGC 4203. In Section 2 we describe the data (Hi, 
EUV, infrared as well as ancillary optical imaging) and the method¬ 
ology used for our work. In Section 3 we discuss global properties 
of the H I, EUV, infrared and optical images, present the derivation 
of star formation rate (SER) from the data and compare the SE ef- 
flcieny in the Hi disc with that of late-type and dwarf galaxies. In 
Section 4 we discuss the origin of the H i gas and its impact on the 
host galaxy. In Section 5 we give a summary of the work. 
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2 DATA AND ANALYSIS 

2.1 WSRT data reduction 

We select NGC 4203 from the sample of ETGs with extended H i 
discs presented by S12. We observe this galaxy for 9x12 hours us¬ 
ing the Westerbork Synthesis Radio Telescope (WSRT). Compared 
to the original 12-h observation by SI2, this longer integration al¬ 
lows us to image the H i down to a similar column density but at a 
significantly higher angular resolution. This is desirable to compare 
the HI image to the UV image used to study SF (see below). 

We reduce separately the 9 Hi datasets (12-h integration each) 
following standard techniques with the MIRIAD package ( |Sault,| 
[Teuben & Wright] 1995| l. Having calibrated and subtracted the con¬ 
tinuum from the datasets, we combine them in the UV plane and 
image with robust=0.4 weighting, obtaining an Hi cube with an 
angular resolution of 29.3x16.7 arcsec^ (PA=0.5 deg). The veloc¬ 
ity resolution is 16 km s“^ after Hanning smoothing. The noise level 
of the data cube is 0.12 mJy beam"k 

We create an H i image by summing flux at all velocities us¬ 
ing only pixels included in a mask. In order to include diffuse, faint 
emission in the mask, we select pixels above a certain threshold in: 
(i) the original data cube; (ii) four data cubes smoothed spatially 
with a 2D Gaussian of FWHM 30, 45, 60 and 90 arcsec, respec¬ 
tively; (iii) three data cubes smoothed in velocity with a Hanning 
filter of width 32, 64 and 112 kms"^; (iv) six data cubes obtained 
smoothing in velocity the 30- and 60 arcsec resolution cubes as in 
point (iii). We create the mask using a 3cr threshold except for the 
90-arcsec resolution cube, for which we adopt a 4cr threshold. 

We determine a detection limit from the original data cube 
as a 5cr signal in a single velocity resolution element, obtaining 
2.3 X 10^^ cm"^. We obtain the final Hi image by applying the 
WSRT primary beam correction to all pixels above this detection 
limit in the moment zero image derive from the masked H i cube. 
We show the Hi image in Fig. ^ (right panel) and a zoom in on the 
central 300 arcsec in Fig.(top-left panel). 

2.2 GALEX data reduction 

We obtain GALFX UV images of NGC 4203 from the Mikulski 
Archive for Space Telescopes. We use data from the Calibration 
Imaging, whose goals and specifications are described in|Martin &| 
IGALFX Team] ( [2005] ). The GALFX field of view (FOV) is 1.2 
degrees in diameter and the instrument has two main bands: the far- 
UV (FUV) and the near-UV (NUV). The effective wavelengths for 
the FUV and NUV channels are 1516A and 2373A, with 4.5 arcsec 
and 6 arcsec resolution (FWHM), respectively. Details can be found 
in [Morrissey et al.j ( [2005 j [2007] ). Our study is mostly based on the 
FUV data because the FUV emission is more sensitive to recent 
star formation than the NUV emission ( [Bianchi|2011[ ). 

Since NGC 4203 has not been observed as a target galaxy, the 
position of the galaxy is different in all the 10 GALFX images. 
Therefore, we have selected 5 GALFX images where NGC 4203 
is sufficiently far from the edge of the GALFX FOV to grant good 
image quality and reasonable background uniformity. The total ex¬ 
posure time of these 5 UV images is 1186 sec. For the purpose of 
comparing the distribution of UV emission and Hi gas, we anal¬ 
yse a 1500 X 1500 arcsec^ GALFX image centered on NGC 4203. 
Figure shows the central 300 arcsec of the FUV image, which is 
the only region where emission from SF in the H i disc is clearly 
visible by eye. 

We make use of the SDSS photometric catalogue ( |Ahn et ar[ 


| 2012 | ) to mask background and foreground objects based on their 
type (unresolved or extended), brightness, colour and Petrosian ra¬ 
dius. For the unresolved sources we adopt a radius equal to the 
GALFX point spread function (PSF), otherwise we use the Pet¬ 
rosian radius of the object itself. Following visual inspection of the 
UV images, we choose to mask the objects with g' -r' <1 and 1.5, 
and with g < 19.5 and 21.1 mag for the unresolved and resolved 
sources, respectively. We estimate the local background from the 
masked images by applying a moving mean filter of size 300 x 300 
pixel^ (450 X 450 arcsec^). The filtering is performed iteratively by 
clipping three times at a 4 . 5 cr level. To avoid over-estimating the 
background, we mask a large area (300 arcsec diameter) around 
the centre of the galaxy. We subtract this background from the orig¬ 
inal, masked image. Additionally, a flat residual background is sub¬ 
tracted from the UV image before proceeding with our analysis (see 
Sec.[T3]). 

Finally, we correct the FUV image for the effect of Galac¬ 
tic extinction. We use an F(B-V) = 0.012 based on the dust im¬ 
ages by [Schlegel, Finkbeiner & Davisj ([1998[) and assume that 
Afuv=8.24 X i'CB-V) and Anuv=8.2 x E(B-Y) ( |Wyder et al.|2007^ 
Thus, we use the following equations for the Galactic correction: 

huV,corr = huV X 1.1 (1) 

^NUV,corr = ^NUY X 1.08 (2) 

2.3 Ancillary optical and infrared images 

We support our study of SF, which is primarily based on the UV 
images described above, with ancillary images at 8 jum, 24 fim and 
optical wavelengths. The 8 fim mid-infrared emission is generally 
considered to trace the PAHs heated by UV photons ( [Leger & Puget| 
|1984|[Seilgren|1984[ ). It is thought that 2/3 of the 8 //m emission 
from a galaxy originates from heating by stellar populations 100 
Myr or younger ( jCalzetti et al.|2005[[Calzetti|2013| ), and that PAH 
emission is stronger for increasing metallicity and dust content of 
the interstellar medium (e.g., [Fngelbracht et al. [2005[[Marble et al.| 
|2010[[Calzetti|2013] l. The 24 fim emission is generally produced by 
dust following absorption of UV photons emitted by young stars, 
and traces SF on a time scale of 10 Myr (jCalzetti et al. [2005 1 [Leroy | 
[et al.|2008| ). 

In this paper we use archival 3.6 fim, 8 jum and 24jj.m images 
from [Fazio et ar[ ( [2004| ) obtained with the Infrared Array Cam¬ 
era and The Multiband Imaging Photometer on the Spitzer Space 
Telescope ( [Rieke et al.|2004|[Wemer et al.|2004| ). We remove the 
background from each image by subtracting the average pixel value 
calculated in a number of empty regions. 

The 8 yum image includes a contribution from stars in the 
galaxy, which we remove by subtracting the 3.6 pm image scaled 
by a factor X 3.6 as in [Shapiro et al.[ ( |^10| ). These authors find that 
X 3 6 has very low scatter about a median value of 0.264 for FTGs, 
and we adopt this value for NGC 4203. The bottom-left panel of 
Fig. shows the resulting 8 pm-non stellar image. Similarly, the 
24 pm flux is also contaminated by emission from old stellar pop¬ 
ulations and hot, circumstellar dust (e.g., [Salim et al.|200^[Kelson| 
|& Holden[|2010[ [Martini, Dicken & Storchi-Bergmann|[2013| ). To 

remove this emission, we apply the same scaling method used for 
the 8 pm. We have adopted a scaling factor from [Martini, Dickenj 
[& Storchi-Bergmann[ ( [2013[ ) that study the empirical 24 pm scaling 
factor for dust-free FTGs. 

We also use deep optical images of NGC 4203 that are ob¬ 
tained in the g' and r' bands using the MegaCam camera on the 
Canada-France-Hawaii Telescope. These observations are obtained 
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Figure 2. Zoomed-in images of the deep H i (top-left), GALEX-FUV (top-right), Spitzer Sjum infrared (bottom-left) and deep MegaCam g' - r' colour (bottom- 
right) images. All images show a field of size 300" x 300" (~ 22kpcx22kpc). Dark areas correspond to redder regions in the deep MegaCam colour image. The 
red contours on the Hi and FUV images represent an Hi column density of 1.8 x 10^^ cm“^. The location of the inner Hi ring is shown by the green-dashed 
lines. 


as part of the MATLAS / ATLAS^° Large Program. The observing 
strategy and data reduetion procedures with the LSB-Elixir pipeline 
optimized for the detection of extended low surface brightness fea¬ 
tures are described in |Duc et "^ ( 120151 ). Further processing is car¬ 
ried out to remove the extended artificial halos surrounding two 
nearby bright stars. We determine the stellar PSFs, including their 
extended wings, with a manual procedure (as explained in |Duc| 
|et aL|2015] ). On the resulting image (see Fig.[^ top), the early type 
galaxy exhibits a rather regular extended outer stellar halo, and a 
dwarf galaxy can be seen east of NGC 4203. 


We compute a model of the NGC 4203 with the ellipse fitting 
procedure of IRAF for both the g' and r' band images and then sub¬ 
tract these models from the original images. The resulting residual 
image together with the g' band image of the dwarf galaxy is shown 
at the bottom panel of Fig. and it is obvious that some extended 
spiral/ring structures show up in the main body of the lenticular 
galaxy. 
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Figure 3. Top: Composite g' + r' Megacam image of NGC 4203 with the dust lanes superimposed in red. The dust lanes are obtained by summing the residual 
g' + r' images after a galaxy model subtraction (see Sec. |2.3) . This image also shows the disrupted companion dwarf galaxy at the east of the galaxy. Bottom: 
MegaCam g' band image of the disrupted companion dwarf galaxy and residual map of NGC 4203. It exhibits a possible faint tail East of our galaxy which as 
the same direction as the two tidal tails of the companion dwarf (see Sec. |3.1) . The white bar shows a lenght of 1 arcmin. 


2.4 SFR calculation 


We estimate the SFR using the various tracers described above 
(FUV, non-stellar 8 jjm and 24 jjm). In order to convert FUV emis¬ 
sion into SFR we use a conversion from |Bigiel et al.| ( [20I0| ), 

2sFR[MGyr-‘kpc“^] = 0.68 x 10“^* x 7Fuv[erg s“‘Hz-‘kpc“^], (3) 


where Ipuv is the FUV intensity per unit area and the initial mass 
function (IMF) is assumed to be of a Kroupa type. This estimate 
could in principle be corrected for internal extinction on the basis 
of the Hi image using t he relation N(ni)IE{B - V) - 5 X 10^^ 
cm"^mag"^ presented by Bohlin, Savage & Drake (19781. This re¬ 


lation is derived for the Milky Way and therefore is likely to provide 
an upper limit on the extinction in metal poor environements. In the 
case of NGC 4203, the peak H i column density is ~ 4.5 x 10^*^ cm"^ 
and, thus, the maximum extinction should be around ~ 25 percent. 
In the outer regions, where the H i column density is much lower 
(~ 3 X lO^^cm"^), the extinction decreases to a negligible ~ 2 per¬ 
cent. In this paper we do not apply this correction but instead trace 
SF in dusty environments by using infrared imaging. 

Although we do not expect a large amount of dust in the outer 
regions of NGC 4203, Fig. shows that in the inner regions dust 
is present and this will affect the calculation of SFR. As mentioned 
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in Section [O] the FUV emission absorbed by dust is re-emitted at 
infrared wavelengths. We use the conversion from |Wu et al.| ( [2005] ) 
to estimate SFR from the 8 jim non-stellar image: 

2^SFR[Moyr"^kpc"^] = 8.43 x 10"^/8;,«(non stellar), (4) 

where is the non-stellar 8 fim intensity in MJy ster"^ 

We also combine the FUV and non-steallar 24 jjm images to 
estimate the SFR as in |Leroy et al.| ( |200^ : 

2^SFR[Moyr"^kpc"^] = 8.1xlO"^/Fuv+3.2xlO"^/24^;„(non stellar), (5) 

where /puv and are in MJy ster"^ and similar to Eq.j^a Kroupa 
type IMF is used. 

Below we study the SF at the original resolution of the FUV 
and infrared images as well as at the resolution of the H i image in 
order to compare the distribution of neutral gas and SFR. 

3 GLOBAL PROPERTIES OF FUV, H i, AND ANCILLARY 
IMAGES 

3.1 Flux distribution of the images 

Early studies emphasize that the H i properties of NGC 4203 are 
similar to those of spiral galaxies. Eor example, [Bur stein & Krumm| 
( |1981| ) and [van Driel et al.| ( |1988| ) find a distorted rotation curve 
from the H i velocity field, possibly suggesting an oval distortion 
or a kinematical warp jBosma |1978| ). They find that the Hi disc is 
3 times larger than the optical body and this is confirmed later by 
S12 using deeper data. 

As described in Sectionj^ our new Hi image has similar sen¬ 
sitivity as that of S12 but a much better resolution. We detect Hi 
down to a column density of 2.3 x 10^^ cm“^ at a resolution of 29.3 
X 16.7 arcsec^. Using the new H i observations, one can easily see a 
dense gas ring whose size is similar to that of the stellar body (Eig. 
[^. This ring is embedded in a much larger Hi disc characterised 
by long spiral-like arms. The H i disc extends out to ~ 40 kpc from 
the galaxy centre and contains ~ 9 times more H i gas than the inner 
ring. Tablej^lists some basic galaxy properties derived from the H i 
image. 

3.2 Impact of the H i on the host galaxy 

Eigurej^shows that NGC 4203 hosts some diffuse EUV emission at 
the location of the inner H i ring. Additionally, at the same location, 
we see spiral-like structure in both the 8 jim and g' - r' images, 
which is clumpy and almost identical in these images. Moreover, 
the extended spiral/ring structures in the model subtracted deep op¬ 
tical image (bottom panel of Eig. correspond to the dust lanes 
visible in the MegaCam g' - r' colour image. 

In contrast, the bright H i arm ~ 2 arcmin west of the galaxy 
(visible in the same figure and referred to as western Hi arm in 
the rest of this paper) does not exhibit any clear EUV emission. A 
similar picture is evident from the 8 fim-non stellar image and the 
g' - r' image. While there is clear emission from PAHs and dust 
absorption in the H i ring, these features are absent from the region 
occupied by the western Hi arm. The question is why these two 
regions have such different levels of the EUV and 8 /rm emission 
despite the very similar H i column density. We will come back to 
this point in Sec |4.1| 

As indicated in Sec. |2.3| the deep optical images show a dwarf 
satellite ~ 7.2 arcmin (32 kpc) east of the NGC 4203. This dwarf 
is red and at a similar redshift with Vheiio - 1067 + 24 kms"^ 


Table 2. Observational results for NGC 4203 


Quantity 

Value 

Unit 

_l_b,c,d 

logM(Hi) 

9.12 

Mo 

2% 

\0gM(ni)Ring 

8.12 

Mo 

5 % 

\0gM(lii)Disc 

9.07 

Mo 

2% 

\0g[M(ni)/M,tellar] 

-1.442 

- 

- 

FUV^ 

16.10 

ABmag 

0.05 

FUV(R<100”) 

16.92 

ABmag 

0.06 

SFR^/„g^8jUOT 

8.4x10-3 

Mo yr"i 

3.0x10-3 

^^^Ring,FUV+24fim 

4.3x10-3 

Mo yr"^ 

2.0x10-3 

^^^Ring,FUV 

3.8x10-3 

Mo yr"i 

1.0x10-3 

SF^Disc 

9.0x10-'^ 

Mo yr"^ 

4.0x10-'* 

S¥^Rifig+Disc,FUV 

4.7x10-3 

Mo yr"^ 

1.0x10-3 

Zy FRRing+Disc,FUV 

5.9x10-3 

Mq yr"^kpc"^ 

1.3x10-3 


Note: 

^ We take into account the Galactic extinction correction for all the values 
below. Look at the text for the details. 

^ The errors are in the same unit as the values. 

The errors in Hi mass values are derived from the noise level (Icr) of the 
mask file that is used to create the H i image. We use an annulus that has the 
same size of the H i beam. 

^ The errors in the SFR and FUV magnitude values are calculated from the 
uncertainties in the photometric measurements of the FUV, 8 //m and 24 urn 
images (Calzetti et al.|2007||Morrissey et al.|2007][Leroy et al.|2008) . 

( jAdelman-McCarthy et al.|2007| ). The stellar mass of the dwarf is 
roughly ~ 10^ M© from the SDSS-/ band photometry. 

The bottom panel of Fig. [^reveals prominent tidal tails asso¬ 
ciated to the dwarf and a faint stream, which has the same direction 
as these tidal tails. From this point, it is clear that the satellite dwarf 
galaxy is interacting with NGC 4203, for example, these features 
may be associated with a single on-going minor merger event. We 
will discuss this point in Sec. |4.2| 


3.3 H I, FUV and infrared radial profiles 

We derive the radial profile of H i, UV and infrared emission along 
circular annuli whose width increases logarithmically from inside 
out. The inner region of the radial profiles (1 Reff) is contaminated 
by old stars in the bulge and, therefore, we do not include it in our 
SFR calculation. 

Figure [^ shows the resulting H i and UV radial profiles to¬ 
gether. The first peak of the H i profile indicates the location of the 
inner H i ring ~ 30-80 arc sec from the centre. We find that the slope 
of the declining FUV profile becomes flatter in this region con¬ 
sistent with the presence of diffuse FUV emission associated with 
the star formation (see Fig.[^. Additionally, the FUV-NUV colour 
becomes significantly bluer at the same location. This might be be¬ 
cause the young stellar populations are dominant in this region. The 
second peak in the H i profile of Fig.|^(~ 120-150 arcsec) is caused 
by the high column density gas in the western Hi arm (see Fig.[^. 
Unlike for the first peak of the Hi profile, here the FUV profile 
does not seem to depart from the general declining trend. However, 
the FUV-NUV colour seems to become bluer at this location. The 
most significant diffuse FUV emission in NGC 4203 comes from 
a region close to the UV-bright galaxy bulge. The emission from 
the bulge may contaminate a large region around it because the 
GALEX PSE is characterized by broad wings out to at least 60 arc- 
sec. To quantify this effect we create a bulge model with a Sersic 
surface brightness profile and convolve it with the EUV PSE. We 
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Figure 4. Comparison of different SFR radial profiles together with the Hi profile. We calculate the SFR with Eqs.[^|^and[^ Although, we show the central 
emission from all the bands, we do not include these emissions in our SFR calculation. The dashed part of the profiles traces emission from the bulge and this 
is not included the calculation of SFR in the rest of paper. The gray shaded area and H i profile are the same as in Fig.[^ 


find that the average contamination from the bulge is ~ 29.2 AB- 
mag arcsec"^ between 50 and 60 arcsec. This contribution to the 
total emission is very low and we ignore this effect in our SFR cal¬ 
culation. 

In order to take into account the significant 8 jjm non-stellar 
emission in the inner H i ring (see Fig.[^, we derive the SFR radial 
profiles based on this 8 fim emission as well as the 24 fim emission 
(see Sec. |2.3| and Eqs. |^and|^. We show the resulting SFR radial 
profiles together with the FUV one in Fig. All profiles show a 
similar trend as their slope changes at the location of the H i ring 
(in fact, at 8 jjm the profile increases significantly between 30 and 
60 arcsec). The infrared based SFR profiles are up to a factor ~ 3 
higher than the profile based on the FUV emission alone. This dif¬ 
ference is larger than the errors in the photometric measurements 
and is likely due to infrared emission from obscured star formation. 
The photometric uncertainty for the Spitzer 8 and 24 pm images is 
~ a factor of 2, and this uncertainty is also similar for the FUV 
measurements ( |Calzetti et al.|2007||Salim et al.|2007| l. However, at 
present, the depth of the Spitzer images and the small FOV imply 
that the UV images are the only reliable tracer of SF at larger ra¬ 
dius. Table [2] shows the SFR efficiencies at the location of the inner 
Hi ring for all 3 tracers as well as their approximate uncertainties. 

3.4 SFR in terms of H i column density 

In Fig.j^we compare H i and SFR surface density by plotting the H i 
and FUV radial profiles of Fig. [^against one another. In this figure 
points are colour-coded according to the distance from the galaxy 
centre. The reason that the track described by our measurements is 
not monotonic is the result of the two H i peaks (one corresponding 


to the inner Hi ring and the other to the outer Hi arm; see Figs.|^ 
an while the UV radial profile decreases monotonically with ra¬ 
dius. For comparison, the yellow area shows the spatially resolved 
relation between H i and SFR surface density inside r 25 for a sam¬ 
ple of nearby (D < 12 Mpc) spirals covering a range of morpholo¬ 
gies ( [Bigiel et al.|2008| ). In this inner region, interstellar medium of 
these galaxies is dominated by molecular hydrogen. The green ver¬ 
tically and horizontally shaded areas show the same relation for the 
outer parts (i.e., 1-2 x r 25 ) of dwarf and spiral galaxies, respectively 
( |Bigiel et al.|2010| ). Details about their sample, e.g., the THINGS 
sample, are given in [Walter et al.| ( [2008] ). This figure shows that 
the SF efficiency in the bright H i ring around NGC 4203 is consis¬ 
tent with that of late-type galaxies. Moreover, the SF efficiency in 
the outer regions of NGC 4203 is also consistent with that of outer 
regions of spiral and dwarf galaxies. 

So far we have compared H i and SFR surface density at fixed 
radius. Below, we perform a pixel by pixel comparison as imple¬ 
mented by previous authors (e.g., [Bigiel et al.|2008||2010l|Boissier| 
[et al.[201^ . To do so, we smooth and re-grid the FUV and infrared 
images to the resolution and coordinate grid of the H i image. We 
find this comparison to be characterised by a very large scatter. For 
this reason, we bin the points along the HI column density axis, 
and show the average SFR surface density of each bin as well as 
the scatter in Fig.[^ 

This plot shows a weak trend between Hi and SFR surface 
density as well as a large scatter, represented by the blue bars. This 
large variation of SF efficiency at fixed H i column density is also 
clear when one compares the H i ring and the western H i arm men¬ 
tioned above. The green asterisks and red pentagon indicate the 
average values of SFR surface denstiy and Hi column density in 
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Figure 5. Top panel: Radial profiles of the H i (blue) FUV (red) and NUV (yellow) emission. The left y-axis shows the UV magnitude and the right y-axis 
shows the Hi surface density. The green-vertical line shows the effective radius of the galaxy. The gray-shaded region indicates the inner Hi ring. Bottom: 
Radial FUV-NUV colour profile. The error bars represent the uncertainty in determining the background for both images. 


these two regions, respectively, and are calculated on the basis of 
the FUV image. Combining the 24 pm image with the FUV image, 
the SFR estimate increases further, as illustrated by the first green 
line. If we calculate the SFR with the 8 pm non-steallar emission, 
we obtain the value that is indicated by the second green bar. As 
argued above, this strengthens our claim that SF efficiency in the 
inner H i ring is as high as that of spiral galaxies, while it is much 
lower in the western H i arm despite the similar gas density. How¬ 
ever, we note that even if the SFR is low in the western H i arm, this 
SFR is still consistent with the outer regions of spirals and dwarf 
galaxies. 


4 DISCUSSION 

4.1 Star formation efficiency 

The previous section demonstrates that at any given Hi column 
density there is a large range of SFR surface densities and, there¬ 
fore, SF efficiencies in the Hi disc of NGC 4203. This is high¬ 
lighted by the comparison between the H i ring and the western H i 
arm presented in Fig. [ 7 ] There are a number of possible explana¬ 
tions for the difference between these two regions: (i) there may be 
no direct relation between Hi and SF; (ii) the Hi column density 
of the arm might appear to be higher because of projection effects 


(e.g., the Hi arm stretches away from the plane of the disc), and 
the actual density of the gas could be too low to form stars; (iii) the 
temperature of the H i in the arm may be higher than in the ring, 
making it more difficult to form stars; (iv) The metallicity and dust 
content of the inner and outer regions might be different due to a 
difference in formation processes; (v) the western Hi arm may be 
stable against gravitational instabilities, preventing the formation 
of stars. 

Both the second and third points would result in a broaden¬ 
ing of the HI spectral profile at the location of the H i arm. This 
is, however, not observed, leaving us with no clear indication in 
favour of these options. On the other hand, the fourth point seems 
an attractive explanation. [Miller, Knezek & Bregmanj ( |1999| ) ob¬ 
serve a quasar ~ 2 arcmin South-East of the galaxy. Interestingly, 
they detect one of the closest known damped Lycr (DLA) system. 
This DLA system belongs to NGC 4203, and the metallicity of 
the gas to be [Fe/H] = -2.29 + 0.1. They also show a number of 
heavy-element abundances such as [Mn/H]< -0.68 and [Mg/H] >- 
2.4, respectively. The velocity of the lines are identical to that of 
the HI emission from the disc of NGC 4203 at the position of the 
quasar within (+10 kms"^. Although the metal abundance values 
are very low, the depletion of the heavy elements onto small dust 
grains can affect the observed measurements. |Vladilo| ( |2002| ), for 
instance, show that the corrected [Fe/H] abundance of a DLA sys- 
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Figure 6. Comparison between Hi and FUV surface brightness obtained plotting the radial profiles of Fig. [^against one another. The triangles indicate Hi 
detected regions with no detectable FUV emission, we estimate the upper limit as equal to the r.m.s. of the FUV background image. The colour bar shows the 
distance from the centre of the galaxy. The yellow area shows the relation between Hi and SFR in the sample of late-ty pe galaxies of|Bigiel et al.H2008) . The 
vertical and horizontal green-shaded areas show the location of dwarf galaxies and outer regions of spirals, respectively ( jBigiel et al.|2010) . In this Figure, the 
left axis shows the FUV surface brightness and the right axis shows the SFR calculated with Eq.[^ 


tern can be on average 0.5 dex (in some extreme cases 1 dex) higher 
than that of observed. Even if we eonsider 1 dex depletion effeet, 
the metallicity of the outer Hi disc remains very low (i.e., [Fe/H] ~ 
-1.29). 


This is likely quite different from the inner Hi ring, which 
is known to contain a considerable amount of dust (Fig. [^. In 
particular, we know that the PAH moleeules traced by the 8 fim 
emission are vunerable to highly energetic photons in dust poor 
environments due to low dust shielding (e.g., [Boselli, Lequeux &| 
|Gavazzi|2004[ |Calzetti et al.|2005| l, and the PAH emission gener¬ 
ally vanishes in regions with metallicity below solar ( |Engelbracht| 
|et al.|2005) l. Thus the metallicity in the inner region should be sig¬ 
nificantly higher than in the outer regions. Consequently, different 
formation processes for the inner and outer H i regions can explain 
these different metallicities. 


The last possibility is the differences of the gravitational insta¬ 
bilities of the Hi gas in the inner and outer regions. For example, 
star formation can occur in regions where the gas disc is unstable 
against large-scale collapse. If the gas is stable against it, Coriolis 
forces counteract the self-gravity of the gas and suppress star for¬ 
mation ( [Leroy et al.|2008| and references therein). This might ex¬ 
plain the different SF efficiencies. However, in order to understand 
this effect, we need further investigation. 


4.2 Origin of the H i gas 

There are different scenarios that could explain the origin of the H i 
disc of NGC 4203. One possibility is a high-angular-momentum, 
gas-rich major merger. Such a merger could result in the formation 
of a large gas disc (e.g., [Barnes |2002| ). At the same time, it would 
not necessarily trigger the central burst of SF usually thought to oc¬ 
cur during mergers (e.g., [Pi Matteo et al.[2007[[Serra et al.[2008[ l, as 
evidenced by the old stellar populations in the centre ofNGC 4203. 
If this scenario is correct, we should able to see a residual from the 
major merger close to the eentre, or a distortion in the outer stel¬ 
lar halo of NGC 4203 in the deep optical image. However, we do 
not detect any of these features in Fig. Another problem with 
the major merger scenario is the extremely low metallicity of the 
outer disc which may rule out that the H i once belonged to a large 
galaxy. 

The other possibility is that the gas has been accreted from 
the inter-galactic medium. The disturbed appearance of the H i dise 
may be explained by some minor interaction with a satellite. In¬ 
deed, as explained in Sec. O there is a nearby dwarf satellite, 
which is most likely interacting with our galaxy. Although this 
dwarf is positioned at the end of one of the H i arms of NGC 4203 
(Fig. [^, we do not detect any Hi emission at that position and, 
therefore, rule out that it is the donor of the Hi gas. 

However, there is another possibility that this gas can be cap¬ 
tured from another gas-rich dwarf galaxy in a single encounter. 
Since the Hi disc is almost regular, this event would happen at 
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Figure 7. SFR surface density as a function of H i column density based on the pixel by pixel method. Here, the FUV pixels are binned as a function of 
Hi column density. The light blue-shaded area shows the scatter of FUV pixels in the bins. The error bars are similar to those in Fig. I^and represent the 
uncertainty in determining the background. The gray error bar on the bottom-right shows the result of same analysis for all the pixels without H i emission. 
The green asterisk and red pentagon show the inner H i ring and the western H i arm, respectively. The two green lines above the green asterisk represent the 
SFRs in the inner H i ring that are calculated from Eq.|^(for FUV-i-24 jum) and|^(for 8 jum), respectively. The blue points represent the upper limits, calculated 
by using the r.m.s. of the FUV background image. The yellow area and green lines are the same as in Figure]^ 


least a couple Gyr ago, making it difficult to detect. This is con¬ 
sistent with the conclusion of [van Driel et al.| ( [1988| ) based on a 
study of the environment around the galaxy. 

In conclusion, accretion from the inter-galactic medium would 
explain the low metallicity of the outer disc, while the inner Hi 
ring may have been enriched of metals by stellar mass loss in the 
galaxy. Moreover, the disturbed appearance of the outer disc can be 
explained by an interaction with the satellite visible in Fig.^an d|3] 

We have long known that a large fraction of all ETGs host 
low levels of on-going or recent SF. This SF activity is detected in 
their centre using optical spectroscopy (e.g., |Trager et al.|2000| | or 
UV photometry (e.g., |Yi et al.|2005| l. Similar to extended-UV discs 
in late-type galaxies, star formation can also be detected at large 
radius in early-type galaxies using UV and deep optical imaging 


(e.g.. 

Salim & Rich|2010[|Salim et al.|2012[|Eang et al.|2012[|Duc| 

|et al. 

2015|l. In a few cases it has been possible to establish a link 


between the presence of Hi and SF at the outskirts of ETGs (e.g., 
Cortese & Hughes|2009[[Donovan et al.|2009[[Thilker et al.|2010[ 

Due et aL|2014| l. Here we explore such link for NGC 4203, and in 
this section we quantify the impact that the SE hosted within the H i 
disc has on the morphology of the host. 

Table [2] shows the result of our SER calculation based on the 
various tracers discussed above, and performed separately for the 
inner H i ring and the remaining, outer H i disc. In particular, the 
inner H i ring hosts a significant amount of dust-obscured SE, and 


the SER is best estimated using a combination of 24 jim and EUV 
emission. 

The inner H i ring and outer H i disc are very different from 
each other. The former contains only ~10 percent of the total Hi 
mass but hosts ~97 percent of the total SE. The latter hosts the 
remaining ~90 percent of the H i mass but only ~3 percent of the 
total SE. As a result the H i depletion time is much longer in the 
outer disc (~ 10^ Gyr) than in the inner H i ring (~5 Gyr). If nothing 
happens to the galaxy, the H i disc could survive for hundreds of a 
Hubble time at the current SER. 

Overall, the detected amount of SE has a negligible ef¬ 
fect on the stellar disc growth and, therefore, the morphology of 
NGC 4203. Based on the total stellar mass reported by |Cappellari| 
|et al.lpOT^ and on the B/D ratio given in |Krajnovic et al.| ( 20I3| l, 
the mass of the stellar disc is ~2.7 xlO^^ M©. At the current SER, 
this mass would grow by a mere 0.1 percent in 1 Gyr. The star for¬ 
mation rate per unit stellar mass in the stellar disc is of just ~ 10"^^ 
yr"^ 


If we suddenly convert all the gas in the outer Hi disc of 
NGC 4203 to stars, the R band surface brightness of this region 
would be 26.8 mag arcsec“^. This is ~ 2 magnitude lower than Ma- 
lin 1, a giant low surface brightness galaxy that has a low surface 
brightness and very large H i disc ( |Lelli, Eraternali & Sancisi|20T0] |. 
In conclusion it is unlikely that NGC 4203 will become a galaxy 
similar to Malin 1. 


















































Star formation in the outer regions ofNGC 4203 11 


5 CONCLUSIONS 

NGC 4203 has an H i disc that consists of two separate components. 
The first component is the inner star forming H i ring. This ring con¬ 
tains metal-rich gas and a large amount of dust. The SF efficiency in 
the HI ring is comparable to that of inner regions of spiral galaxies 
and the H i depletion time is of the order of a few Gyr. The second 
component is the outer H i disc, where the metallicity is likely to 
be low. The SF efficiency in this disc is much lower than that of 
centre of spirals even when the Hi column density is high (e.g., in 
the western Hi arm) and the Hi depletion time is of lOO’s of the 
Hubble time. However, the SF effieciency in the outer Hi disc is 
still consistent with that of the outer regions of spiral galaxies and 
dwarfs. 

These two components of the H i structure may be formed in 
different ways. On the one hand, the H i ring is most likely enriched 
by stellar proeesses within the galaxy. On the other hand, the outer 
HI disc might be accreted from the low metallicity IGM because 
the outer disc is formed with metal-poor, high column density gas 
and we cannot find a suitable donor. Subsequent interaction with 
satellites, which is possibly on-going, may explain the current dis¬ 
torted morphology of the H i disc. 

We also deteet spiral-like strueture in the 8 fim and deep g' - r' 
images. These spiral-like structure is extending to 3 effective radii 
(~ 6kpe), and, because of its dumpiness, is likely organized by the 
star formation process. 

Finally, despite the large Hi reservoir and the detection of 
some SF, the stellar disc of NGC 4203 is growing at a very low 
rate and the morphology of the galaxy is unlikely to change signif¬ 
icantly in the foreseeable future. 
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